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Abstract 



Quark matter at asymptotically high baryon chemical potential is in a color 
superconducting state characterized by a gap A. We demonstrate that al- 
though present weak-coupling calculations of A are formally correct for /i — > 
oo, the contributions which have to this point been neglected are large enough 
that present results can only be trusted for ^ » ~ 10^ MeV. We make 
this argument by using the gauge dependence of the present calculation as a 
diagnostic tool. It is known that the present calculation yields a gauge invari- 
ant result for /i — > oo; we show, however, that the gauge dependence of this 
result only begins to decrease for n ^ fic, and conclude that the result can 
certainly not be trusted for fi < fic- In an appendix, we set up the calculation 
of the influence of the Meissner effect on the magnitude of the gap. This 
contribution to A is, however, much smaller than the neglected contributions 
whose absence we detect via the resulting gauge dependence. 
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I. INTRODUCTION 



The starting point for a description of matter at high baryon density and low temperature 
is a Fermi sea of quarks. The important degrees of freedom — those whose fluctuations cost 
little free energy — are those involving quarks near the Fermi surface. We know from the 
work of Bardeen, Cooper, and Schrieffer that any attractive interaction between the 
quarks, regardless how weak, makes the Fermi sea unstable to the formation of a condensate 
of Cooper pairs. In QCD, the interaction of two quarks whose colors are antisymmetric 
(the color 3^ channel) is attractive. (The attractiveness of this interaction can be seen from 
single-gluon exchange, as is relevant at short distances, or via counting strings or analyzing 
the instanton induced coupling, as may be relevant at longer distances.) We therefore expect 
that under any circumstance in which cold dense quark matter is present, it will be in a 
color superconducting phase pj-^f. The one caveat is that this conclusion is known to be 
false if the number of colors is Ai'c = oo 0]. Recent work |^,^ indicates that quark matter is 
in a color superconducting phase for A^^^ less than of order thousands, and in this paper we 
only discuss QCD with Nc = 3. 

We now know much about the symmetries and physical properties of color supercon- 
ducting quark matter. The dominant condensate in QCD with two flavors of quarks is in 
the color 3a channel, breaking S'f/(3)coior SU{2), and is a flavor singlet Quarks 
with two of three colors have a gap in this 2SC phase, and five of eight gluons get a mass via 
the Meissner effect. In QCD with three flavors of quarks, the Cooper pairs cannot be flavor 
singlets, and flavor symmetries are necessarily broken. The symmetries of the phase which 
results have been analyzed in Ref. [|lOl, and were in fact first analyzed in a different (zero 
density) context in Ref. |Tl[]. The dominant condensate locks color and flavor symmetries, 
leaving an unbroken global symmetry under simultaneous SU{3) transformations of color, 
left-flavor, and right-flavor. In this CFL phase, all nine quarks have a gap and all eight 
gluons have a mass ||10|- Chiral symmetry is spontaneously broken, as is baryon number, 
and there are consequently nine massless Goldstone bosons |]I0[. Matter in the CFL phase 



is therefore similar in many respects to superfluid hypernuclear matter |10,12-|l4|. The fact 
that color superconducting phases always feature either chiral symmetry breaking (as in the 
CFL phase) or some quarks which remain gapless (as in the 2SC phase) may be understood 
as a consequence of imposing 't Hooft's anomaly matching criterion |T^. The flrst order 



phase transition between the CFL and 2SC phases has been analyzed in detail |13,14,16 



but all that will concern us below is that any flnite strange quark mass is unimportant at 
large enough /i, and quark matter is therefore in the CFL phase at asymptotically large /x. 

Much recent work has resulted in two classes of estimates of the magnitude of A, the 
gap in the density of quasiparticle states in the superconducting phase. The flrst class of 
estimates are done within the context of models whose parameters are chosen to give reason- 
able vacuum physics. Examples include analyses in which the interaction between quarks 
is replaced simply by four-fermion interactions with the quantum numbers of the instanton 
interaction [^^|T7| or of single-gluon exchange ||TD|,|T3[ and more sophisticated analyses done 
using instanton liquid models Renormahzation group analyses have also been used 

to explore the space of all possible four-fermion interactions allowed by the symmetries of 
QCD pO| , pT| . These methods yield results which are in qualitative agreement: the gaps 
range from several tens of MeV up to as much as about 100 MeV and the corresponding 
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critical temperatures, above which the superconducting condensates vanish, can be as large 
as about 50 MeV. 

The second class of estimates uses /x — > oo physics as a guide. At asymptotically large /i, 
models with short range interactions are bound to fail, because the dominant interaction is 
due to the long-range magnetic interaction coming from single-gluon exchange [p2| , p3[ . The 



collinear infrared divergence in small-angle scattering via single-gluon exchange results in 
a gap which is parametrically larger at fi oo than it would be for any point-like four- 



fermion interaction 0. Son showed that this collinear divergence is regulated by Landau 
damping (dynamical screening) and that as a consequence, the parametric dependence of 
the gap in the limit in which the QCD coupling g is 

A 1 / 37i'\ 
which is more easily seen as an expansion in g when rewritten as 




5\ng + f{g). (2) 

V ^ 9 

This equation should be viewed as a definition of f{g), which will include a term which is 
constant for g ^ and terms which vanish for g ^ 0. The result (|l]) has been confirmed 
using a variety of methods pi|-Pm , and several estimates of lim^^o fid) exist in the literature. 



For example, Schaefer and Wilczek find p4| , p0 
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(3) 



in the CFL phase (see also Ref. ||2^), and Brown, Liu, and Ren find a result for lim^^o fig) 
which is smaller by (vr^ + 4)/8 — ln2 = 1.04. If this asymptotic expression is applied by 
taking g = g{fi) from the perturbative QCD /5-function (with Aqcd = 200 MeV), evaluating 
A at yU ~ 500 MeV yields gaps in rough agreement with the estimates based on zero-density 
phenomenology. 

The central purpose of this paper is to demonstrate that this nice agreement must at 
present be seen as coincidental, because present estimates for / are demonstrably uncon- 
trolled for g > gc ^ 0.8, corresponding to fi < fic with fi^ ~ 10^ or higher. 

The weak-coupling calculations are derived from analyses (done using varying approx- 
imations) of the one-loop Schwinger-Dyson equation without vertex correction, and (with 
one exception) yield gauge dependent results. However, Schaefer and Wilczek argue that 
the result for lim^^o fio) in such a calculation is gauge invariant. The one calculation which 
is gauge invariant throughout is the calculation of Tc (and hence A since the BCS relation 
Tc = 0.57A holds p5|) done by Brown, Liu, and Ren [^. As in other calculations, however, 
these authors neglect vertex corrections. Our purpose is to use the fact that our calculation 
(like most) is gauge dependent, and only gauge invariant for g —>■ 0, to estimate the g above 
which vertex corrections, left out of all calculations, cannot be neglected. 

We begin by sketching the derivation of the one- loop Schwinger-Dyson equation for A, 
making as few approximations as we can. We solve the resulting gap equation numerically 
in several different gauges. Our results are (yet one more) confirmation of ([l|). Furthermore, 
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we do find evidence that the gauge dependence of / decreases for g — 0. However, this 
decrease only begins to set in for g < 0.8. This implies that the contributions to A which 
have been neglected — like those arising from vertex corrections — only become subleading 
for g ^ gc ^ 0.8. If we translate gc to fic by assuming g should be taken as g{fi), this 
corresponds to fic ~ 10^ MeV. Recent work [^] shows that g should be evaluated at a 
much lower (^f-dependent) scale than fi. This means that the condition g < g^. ^ 0.8 would 
translate into fi > with fic orders of magnitude larger than /Xc ~ 10^ MeV. 

The original purpose of our investigation was to do a self-consistent calculation of the 
influence of the Meissner effect on the magnitude of the gap in the CFL phase. In the 
presence of a condensate, the gluon propagator is modified: some gluons get a mass. In the 
CFL phase, a// gluons get a mass, and this makes a calculation based on perturbative single- 
gluon exchange a self-consistent and complete description of the physics at asymptotically 
large /i, with no remaining infrared problems. (In the 2SC phase, in contrast, the calculation 
of A leaves unanswered any questions about the non-Abelian infrared physics of the three 
gluons left unscreened by the condensate.) We felt that this motivation warranted a self- 
consistent calculation in which we calculate the gap using a Schwinger-Dyson equation in 
which the gluon propagator is modified not only by the presence of the Fermi sea (Debye 
mass. Landau damping) but is also affected by the condensate (the Meissner effect). We set 
this calculation up in an appendix. Previous work, beginning with that of Ref . |Q , shows 
that the form of Eq. (|lD is unmodified by including the Meissner effect, but f{g) is affected. 
Our preliminary results suggest that the changes in f{g) are small, as anticipated in Refs. 
p3| , |2^p7| , |29|j3^ , |33[| . Indeed, the effects of physics left out of the present analysis, which we 
have diagnosed via the gauge dependence of f{g), are much larger than those introduced by 
the Meissner effect at any g we have investigated. 



II. DERIVING THE GAP EQUATION 

In this section, we derive the gap equation for QCD with three massless flavors which is 
valid at asymptotically high densities. We follow Ref. but make fewer approximations. 



Because our point is to stress the importance of effects which we do not calculate, we will 
make our assumptions and approximations very clear as we proceed. In other words, since 
the lesson we learn from our results is that they cannot yet be trusted, it is important to 
detail carefully all points at which we leave something out. 

We use the standard Nambu-Gorkov formalism by defining an eight-component field 
\E' = ip^). In this basis, the inverse quark propagator takes the form 



S-\k)=(^^P' ,T (4) 




where A = 7oA^7o. The color, flavor, and Dirac indices are suppressed in the above expres- 
sion. The diagonal blocks correspond to ordinary propagation and the off-diagonal blocks 
reflect the possibility for "anomalous propagation" in the presence of a diquark condensate. 
We make the following ansatz for the form of the gap matrix P,|lO|,p3,p4[ : 



Ag'(fc) = {XfY\\f).,C^, {At{h)P^{k) + At{h)P.{k)] 

+ (A5)'^^(A5).,C75 {Afiko)PM + Af (fco)P-(^)) (5) 
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Here, a,b = 1,2,3 are color indices, i,j = 1,2,3 are flavor indices, Xf are antisymmetric 
U{3) color or flavor matrices with J = 1, 2, 3, and Aj are symmetric U{3) color or flavor with 
J = 1, . . . , 6, and the projection operators P± are deflned as 

P_ik) = (6) 

with a = 7o7. 

By making this ansatz, we are making several assumptions: 

• First, we have taken A^, A^, Af, and Af to be functions of fco only. All are in 
principle functions of both and \k\, but we assume that they are dominated by 
I A; I ~ fi. This is a standard assumption, and although we do not expect that relaxing 
this assumption would resolve the problems which we diagnose below, this does belong 
on the list of potential cures. 

• Second, we have explicitly separated the gaps which are antisymmetric 3a in color and 
flavor from those which are symmetric 65 in color and flavor and, in both cases, we 
have assumed that the favored channel is the one in which color and flavor rotations 
are locked. The color and flavor structure of our ansatz is thus precisely that flrst 
explored in Ref. |T0[, which allows quarks of all three colors and all three flavors to 



pair. Subsequent work []30| , p2| , |33| , p5| conflrms that this is the favored condensate, and 
we will not attempt to further generalize it here. 

Third, we have assumed that the Cooper pairs in the condensate have zero spin and 
orbital angular momentum. This seems a safe assumption in the CFL phase, where 
the dominant condensate, made of Cooper pairs with zero spin and orbital angular 
momentum, leaves no quarks ungapped. 

Fourth, we neglect ipijj condensates. Since chiral symmetry is broken in the CFL 
phase, these must be nonzero [0]. This applies to both color singlet and color octet 



i/jiIj condensates iSH]. Such condensates are small |T^,|3D|, however, and we expect that 



neglecting them results in only a very small error in the magnitude of the dominant 
diquark condensate. 

The most important assumption we make is that we obtain the gap by solving the one- 
loop Schwinger-Dyson equation of the form 

S-\k) - S^\k) = ^g'J -^TlS{q)TlD:i{k - q) , (7) 

using a medium- mo difled gluon propagator described below and unmodifled vertices 

T.a_f7.AV2 \ 
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Here, 5*0 is the bare fermion propagator with A = 0. Note that we use a Minkowski metric 
unless stated otherwise. We will demonstrate that our results are completely uncontrolled for 
9 > 9c ^ 0.8. This breakdown could in principle reflect a failure of any of our assumptions. 
We expect, however, that it arises because contributions which have been truncated in 
writing are large for g > Qc- That is, we expect that this truncation (and not any of the 
simplifications introduced by our ansatz for A) is the most significant assumption we are 
making. 

We obtain four coupled gap equations 



Afo(A;o) 



-6' 



'6' 



J ^Tr [P±(A:)7, (P+(g)a+(g) + P-iq)a-{q)) 7.] D>^-^ik - q) 



(27r 



j ^Tr [P±(A:)7, (P+(g)6+(g) + P-{q)h.{q)) 7.] D^^{k - q) (9) 



(27r 



where P± means P+ in the Ai equation and P_ in the A2 equation and where 
-Af (go) - A^(go) 



a+(g) 



go'-(l9l+/i)^-4[A^(go) + 2Af(go)]' 

A^(go) - Af (go)l \-ql + {\q\+ fif + (5A^(go) + 7Af (go))(A^(go) + 3Af (go)) 



a_(g) 



-m+ - (A^(go) - Af (go))2] [gg - (|gl + /i)^ - 4(A^(go) + 2Af (go))^] 
-Af(go)-Ai^(go) 



%'-(k1-/^P-4[Af(go) + 2Af(go)]^ 

Atiqo) - Af (go)l [-go^ + (|g1 - /i)^ + {5At{qo) + 7Af (go))(A^i(go) + 3Af (go)) 



b+{q) 



2 - ^|g-1 - /i)2 - (Atiqo) - Af (go))2] [gg - dg] - ^)2 _ 4(A^i(go) + 2Af (go))^] 
Af(go) 



qi-{\q\+f^r-A[A^{qo) + 2Al{qo)f 

A^(go) - Af (go)l [A^(go) + Af (go)! [A^(go) + 5Af (go) 



6_(g) 



b'o -m+ - (A^(go) - Af (go))2] [gg - {\q\ + - 4(A^(go) + 2Af (go))^] 
Af(go) 



%'-(k1-/^)^-4[Af(go) + 2Af(go)]^ 

Atiqo) - Af (go)l lAtiqo) + Af (go)! lA^^qo) + 5Af (go) 



[qi f^)' - (A^^(go) - Af (go))2] [gg - {\q\ - - AiA^qo) + 2Af (go))^] ' 

In a general covariant gauge, the resummed gluon propagator is given by 



(10) 



D^.u{q) 



P 



+ 



p. 



q,iqu 



g2 - G{q) g2 - F{q) q 



(11) 



where G{q) and P(g) are functions of go and |g| and the projectors P^ are defined as 



follows: 



p^ 



'^^^ g2 



(12) 
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The functions F and G describe the effects of the medium on the gluon propagator. If we 
neglect the Meissner effect (that is, if we neglect the modification of F{q) and G{q) due 
to the gap A in the fermion propagator) then F{q) describes Thomas- Fermi screening and 
G{q) describes Landau damping and they are given in the HDL approximation by 




with Qo{x) = - log 



k1 ) \q\ 



x + 1 

X — 1 



(13) 



where = 3g'^fi^/27T'^ is the Debye screening mass for Nf = 3. We discuss the modifications 
of F{q) and G{q) due to the Meissner effect in an Appendix. 

In order to obtain the final form of the gap equation, we need the following trace: 

Tr [P±ik)^, {P+{q)a+{q) + P-{q)a-{q)) 7.] D^^^{k - q) 



r, -lT{k-q)-k{k-q)-q 



+ 



1 =F A; ■ g 



(k-q)l + (k-q)'' 



{k-q)^-F{k-q) 



±2{k~q) ■k{k-q)-q 



(14) 



+a^{q) 



r) —l±{k—q)-k{k—q)-q 
^ {k-q)2-G(k-q) 



+ g^(l±k-q 



, {k-q)l + {k~q)^ 
{k~qr 



{k-q)2-F{k-q) 



=F 2{k — q) ■ k{k — q) ■ q 



(k^qy 



This allows us to recast Eq. (R) into the following form: 



d^q 
(2^)4 



«+WJ M (fc_5)2_G(fc_,) + 



+ j^(l-k-q 



{k-q)^-F{k-q) 

itzS&±^ + 2{k'^q) ■ kik'^q) ■ gM! 



{k-qf 



{k-qf 



. . . , . f, Ak-q)lHk-q)'^ ^lir~\ Uir~\ -i^Z^lo 



{k-q)'^-F{k~q) 

ik-q)lHk^qy q) . k^k^q) . qihsf 



(k-q)^ 



{k-qf 



A^iko) 



(2^)4 



a+(g) 2- 



-l+{k-q)-k{k-q)-q 



{k-qf-G{k-q) 



Vk.q'JlZ,^^-26^yHl^q).q^ 



{k-q)l+ik-g)^ 



{k-qf~F(k~q) 



-2{k — q) ■ k{k — q) ■ q 



+a-{q)\2- 



-l — {k—q)-k{k—q)-q 



(°_,)^ +2(fc-g).fc(fc~g).g^^ 



(fc-g)2-G(fc-f,) 



(A:-g)2-F(fc-g) 

+ 2(A:'^g) ■ fc(A:'^g) ■ q^^ 



(15) 
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Af(A;o) 



6^ /(27r)4 



l — {k — q)-k{k — q)-q 



G{k-q) 



-l-fe-g (t^ +2{k~qyk{k-qyqj^ 



{k-qY-F(k-q) 



+ 




-l+{k^q)-k{k^q)-q , "l+^-g 



(fc-q)2-G{fc-g) 



'2(k-q)-k{k-q)-q 



Af(A;o) 



(2^)4 



(fc-g) 



2 



(fe-q)2 



-l+(fc-g)-fc(fc-g)-g 



-^+k-q TT—T^ 2(k-q)-k(k-q)-q 



+ 



(fc-g)2-G(fc-g) ' 



(fc-g)2 



(fc-g)2 



(fc-g)2-F(fc-g) 



-2(A; — g) ■ k{k — q) ■ q 



(k-qY 



(k-q)- 



-1-k-q 



+2{k-q)-k{k-q)-q 



Ak-q) 







(fc-g)2-F(fc-g) 



(fc-g)2-G(fc-g) 

g)g + 
(fc-g)2 



(15) 



III. SOLVING THE GAP EQUATION 

In order to obtain a tractable numerical problem, we make two further simplifying as- 
sumptions: 

• First, at weak coupling we expect the physics to be dominated by particles and holes 
near the Fermi surface. This manifests itself in Eq. ([I5|) in the fact that a_ and 
h- have singularities on the Fermi surface while a+ and 6+ are regular there, and we 
therefore expect that at weak coupling we can neglect a+ and 6+. Upon doing this, 
we have equations for Aj^ ' which do not involve ' . We are only interested in 

A. S AS 

A^ ' , since Ag ' describe the propagation of antiparticles far from the Fermi surface. 
If we assume that we are at weak enough coupling that a+ and 6+ can be neglected 
(that is if we assume that ''^ -C /i) then we can ignore A^''^ in our calculation of 
Af'"^. (Note that we are not assuming that A^''^ is any smaller than Af'"^; there is 
no reason for this to be true.) We will see that our results break down for g ^ 0.8, at 
which A < 10~^/i. Because A <^ /i, neglecting the effects of A^''^ on Af'"^ should be 
a good approximation, and we do not expect that including these effects would cure 
the problems we discover. This should, however, be investigated further. 

• Second, we set Af = 0, and solve an equation for A^ alone. This assumption is 
in fact inconsistent, as the gap in the symmetric channel must be nonzero. This 
is clear from explicit examination of the gap equations Eq. ([T5|) (and indeed of the 
gap equations of Ref. |]Tn[)- In fact, this result is manifest on symmetry grounds 
[0,0: in the presence of A^ 7^ 0, a nonzero Af breaks no new global symmetries 
and there is therefore no symmetry to keep it zero. Because single-gluon exchange 
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is repulsive in the symmetric channel, this condensate can only exist in the presence 
of condensation in the antisymmetric channel. Explicit calculation [10,30,^| shows 
that the symmetric condensates are much smaller than those in the antisymmetric 
channels. We are therefore confident that keeping Af would yield only a very small 
correction to Z^^. 

We must now solve a single gap equation for A^(/co), which henceforth we denote simply 
as A(A;o). The reader will see below that this equation is still rather involved. Most authors 
have made further approximations, valid for (7 — > 0. Because we make no further approxima- 
tions, our results cannot be gauge invariant. This allows us to test the claim that the results 
become gauge invariant in the limit (7 — 0, and to use the rapidity of the disappearance 
of gauge dependence as this limit is approached to evaluate at what g the contributions we 
have truncated can legitimately be ignored. 

In order to obtain numerical solutions, it is convenient to do a Wick rotation go iQo 
to Euclidean space, yielding the gap equation 



A(*„) = 



(2^ 



A(go) 



-/i)2 + 4A2(go) 



/i)2 + 5A2(go)) 



^(9g + (k1-/i)^ + A2(go))(gg + (|g1- 
1 — {k — q) ■ k{k — q) ■ q 
r {k - q)l + {k- + G{k, -q,^\k-q\) 



/i)2 + 4A2(go)) 



l + k-q 



+ 



+ 2{k -q) ■k{k-q)-q 



(k-q)l 



(k 

-l + k-q 



- q)l + {k- q)' 

-(k-q)l + (k-q)-^ 
(k-q)l + (k-q)^ 



F{ko - go, \k - q\) 
— 2(k — g) ■ k(k — q) ■ q 



{k-q)l + (k-q)^ 



(16) 

{k-q)l + {k-q)^ 

The integral over the azimuthal angle is trivial, and we therefore have three integrals to 
do. We do the remaining angular integral analytically, after making a change of variables. 
We define 



g' 



k 



q 



because the integration over the polar angle 9 is simpler when the momentum integration 
is done over g'. The simplification arises because there is no longer any angular dependence 
in the functions F and G: 

F{ko - go, |A; - gl) = F{ko - go, |g'|) 

and similarly for G. After doing the angular integral, the gap equation reduces to a double 
integral equation with integration variables |g'| (which we henceforth denote g) and go: 



A{ko 



where 



487r3 



dqo / dq 



A(go) 



{ko 
A(go) 



go)^ + g^ + G{ko - go,g) 



iciqQ^q) 



(ko - qoY + g2 + F{kQ - go, g) 



lF{ko,qo,q) +^ 



?A(go 



(ko - go)^ + q' 



;I^{ko,qo,q) 



(17) 
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3?/xV^o+4A2(5o) y^g + 4A^ 

^ 4(gg + A^(go) + q'){q' + V - ql - A^(go)) ^^^^^^ g 

?,qiJ?^ql + A2(go) ^^^Ta^ 
12A^(go) + 6gg-2g^-2V 



lF{ko,qo,q< = 

2((gg + 4A^(go))(fco - gp)^ - q%q^ - V + + 4A^(go)) 



arctan 



^ 4((gg + A^(go))(A;o - gp)^ - - V + gg + A^(go)) ^^^^^^ g 



3g/xVgg + A2(go)((A;o - go)^ + g^) ^g^ + A2(go) 

6g4 + 2(fco - go)^(-2g^ + 12^^ - 3gg - 6A^(go)) 
^l^^{{ko-qor + q') 



Idko,qo,q< ^) = 

2(gg + 4A^(go) - (fcp - go)^)(g^ - + gg + 4A^(go)) 



arctan 



3/iyg2 + 4A2(go)((A;o - go)^ + g^) ^g2 + 4A2(go) 

4(gg + A^(go) - (fcp - go)^)(g^ - 4/.^ + gg + A^(go)) ^^^^^^ g 

3/xygg + A2(go)((A:o - go)^ + g^) ' ^gg + A2(go) 

2g(2g^ - 3(feo - gp)^ - 12^^ + 3gg + 6A^(go)) 
3//2((A;o-go)2 + g2) 

laiqo, q>^J)^ 

{ql + 4A^(go) + g^)(g^ + 4/x^ - gg - 4A^(go)) (^^^^^ q 

2>q^^ ^gg + 4A2(gp) V'"" '''' y'gg + 4A2(gp) 

— arctan . ^ ^ | 
^gg + 4A2(go); 

^ 2(gg + A^(go) + g^)(g^ + 4/.^ - gg - A^(go)) (^^^^^ q 

3g//ygg + A2(go) V'"" ^/^^Ta^ 

g-2/x 



■ arctan 



gg + A2(gp), 



^ 4(gg + A^(go) + g^) gg + A2(go) + 



Sq^i gg + A2(go) + (g - 2^^ 

^ 2(gg + 4A^(go) + g^) gg + 4A2(go) + q' 



3qi, gg + 4A2(go) + (g-2/i)2 

12A2(go) + 6gg - 6g2 - 8//^ - 12//g 
3//g 



lF{ko,qo,q > = 

{{qi + 4A2(go))(fco - go)' - q%q' - V + qI + ^^HQo)) 



cLrcijciri 

— arctan -^^^^=^^^= 
^qi + AA^iqo)) 

2{{ql + A^(go))(fco - q,f - q%q^ - V + gg + A^(go)) / ^^^^^ q 

3qf,^^qi + A^{qo){{h - q,? + q') \ '''' y^^g + A^ 

g - 2/i 



— arctan 



4(g^-(gg + A^(go))(fco-gon g2 + A2(go) + g2 



go' + A2(go). 



3g^(g2 + (A;o - go)^) ql + A2(go) + (g - 2/i)2 

2(g^ - {ql + 4A2(go))(fco - go)') ^o' + ^^\%) + 



3g/x(g2 + [k, - go)2) gg + 4A2(go) + (g - 2/i)2 

6g4 + 2(fco - go)'(6/ig + 4/i2 - Sg^ - 6A2(go)) 



3ijq{{ko - go)^ + q 



•2^ 



I^{ko,qo,q> fi) = 

{ql + 4A2(go) - {k^ - q,f){q' - 4/.2 + g^ + AA\q,)) 



arctan 



S/xVgg + 4A2(go)((A;o - go)^ + g^) V V^o + 4A2(go) 

g-2/i \ 



— arctan 



/g2 + 4A2(go); 

2(gg + A^(go) - {k, - q,f){q^ - 4^^ + gg + A^(go)) / q_ 



Sfi^ql + A2(go)((A;o - go)^ + q') \ ^ql + A^{qo) 

q-2fi 



arctan 



go' + A2(go), 



4(gg + A2(go) - (ko - go)') ql + A'iqo) + q' 



3/i(g2 + {ko - go)2) ql + A^{qo) + (g - 2/x) 

2(gg + 4A2(go) - {ko - go)') gg + 4A2(go) + q' 



2 



3^(g2 + (A;o - go)^) ql + 4A^{qo) + (g - 2/i)^ 

6ql + 12A2(go) - 6(A;o - go)' - 8^^ _ i2^g) 



3yu((A;o - go)^ + 



We have solved the gap equation (p!7| ) numerically for several different values of g and 
several different values of ^. It is convenient to change integration variables from go to In go 
and from g to Ing. We evaluate the g integral over a range gmin < g < 10'' with gmin/A* 
chosen differently for each g in such a way that it is less than 10~^A(0) in all cases. The go 
integral is made even in go (by taking the average of the integrand at go and —go) and then 
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evaluated over a range gomin < qo < 100/i, where we chose gomin = Qmm- We have checked 
that our results are insensitive to the choice of upper and lower cutoffs of the integration 
region. It was probably not necessary to choose and gomin quite as small as we did. It 
is, however, quite important to extend the upper limit of the go and g integrals to well above 
H in order to avoid sensitivity to the ultraviolet cutoff.Q We use an iterative method, in 
which an initial guess for A(/co) is used on the right-hand side of (0), the integrals are done 
yielding a new A{ko), which is in turn used on the right-hand side. The solution converges 
well after about ten iterations. All results we show were iterated at least fifteen times. 

Our results are shown in Fig. ^ Note that the output of our calculation is a plot of 
A(go)//i as a function of go/yU for some choice of g and ^. The only way in which fi enters 
the calculation is to set the units of energy. The values of /i shown in Fig. |l| corresponding 
to each value of g do not come from the calculation. They are obtained by assuming that the 
running coupling g should be evaluated at the scale /i and using the one-loop beta function 
with Aqcd = 200 MeV. We include these values of /i to make comparison with the results 
of Refs. p^j33[| easier. If, as seems quite reasonable, g should in fact be evaluated at a 
5f-dependent scale which is lower than fi, then the values of g at which we have done our 
calculations correspond to larger values of n than shown in Fig. |l] [^. Evans, Hormuzdiar, 



Hsu, and Schwetz have obtained numerical solutions to simplified gap equations describing 
the gap in the CFL phase [^. Their results agree reasonably well with the results of our 
calculation done in ^ = gauge but disagree qualitatively with ours in any other gauge. 



Simply setting ^ = 0, as in Ref. p^j33[] , is not a valid approximation at the values of g at 
which we (and these authors) work. 

How should one interpret the results of a gauge dependent calculation, given that at 
any fixed g one can obtain any result one likes if one is willing to explore gauge parameters 
— oo < ^ < C)0? In the present circumstance, the idea is that we expect this calculation to 
give a gauge invariant result in the g ^ limit. More precisely, if we define 

"A(0) 



fig) = In 



371-2 

+ -^ + 5lng (18) 



then we expect / to go to a ^-independent constant in the g ^ limit. In Fig. |^, we plot 
f{g) in five different gauges. From this figure we learn: 

• For any ^, f{g) is a reasonably slowly varying function of g. This confirms Son's result 
(P and justifies an analysis in terms of f{g). 

• It does appear that limg^o f{g) is a .^-independent constant, perhaps not far from 
the estimate of Ref. [Q, namely limg^o/(fi') = 8.88, or that of Ref. |2^, namely 



lim,^o/(^/) = 7.84. 

If we do a calculation in some fixed gauge, we expect that at small enough g this calcu- 
lation yields a good estimate of the true gauge invariant result. By doing calculations 



The one exception, in which we do find some sensitivity to one of our limits of integration, is 
at g = 3.5576. With g this large, we should perhaps have extended the upper cutoff of the go 
integration to 1000 fi, as the results shown in Fig. 1 below make clear. 
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g = 3.5576 = 400 MeV) x 1 0"^ g = 1 .01 49 (n = 1 0^ MeV) 




FIG. 1. The gap A(go) for five different values of the couphng constant g. In each plot, the 
upper, middle, and lower curves are calculations done using three different gauges ^ = — 1,0, 1. In 
each panel, the range over which the qo integral was done is that shown. 
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FIG. 2. The function f{g), defined in Eq. ([l^), for five different values of the couphng constant 
g. At each g, the points (from top to bottom) correspond to different gauges with = —4, —1, 0, 1, 4 
respectively. Note that the horizontal axis is 1/g and fi increases to the right. At the largest value 
of g, we only show ^ = — 1, 0, 1. In Fig. |l|, we have not shown the A(go) curves for ^ = ±4 because 
in these gauges A (go) is very small or large on the scales of Fig. ||. 

in several gauges, we can bound the regime of applicability of this estimate. We can 
only trust our calculation of f{g) in the regime in which the ^-dependence of / de- 
creases with decreasing g. Our calculation of f{g) is completely meaningless unless g is 
small enough that the curves for different values of ^ are converging. Fig. § shows that 
the gauge dependence of our result for / is about the same for all g ^ 0.8. It is only 
for g < 0.8 that f{g) calculated in different gauges begins to converge. At larger values 
of g our calculation provides no guide whatsoever as to the value of / that would be 
obtained in a complete, gauge invariant calculation including all the physics neglected 
in the present calculation. Even at g = 0.8 the values of A(0) differ by a factor of 
about 400 for gaps with ^ = — 4 and ^ = 4. We could make the gauge dependence 
look even larger by choosing larger values of |^|. Our result does not guarantee that 
the calculation is under control forg < 0.8, but it does guarantee that the result is 
uncontrolled and completely meaningless for g > 0.8. 

IV. CONCLUSION 

We have detailed our assumptions and approximations as we made them. Let us now 
ask which of them should be improved upon if we wish to include those contributions whose 
neglect we have diagnosed via the gauge dependence of our results. Note that g = 0.8 
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corresponds to A//x ~ 10~^. Thus, those contributions to / which we have neglected which 
are controlled when A <^ /i are not responsible for the breakdown of our calculation around 
g ~ 0.8. We believe that the assumptions we made in writing the ansatz (^) and the 
assumptions we made in neglecting Af and A^'"^ all introduce errors which are small when 
A ^ /i. (For example, even though neglecting A2 is a source of gauge dependence, we do not 
expect that remedying this neglect would change f{g) appreciably in any gauge at ~ 0.8, 
where A/n is so small.) Hence, we believe that it is the assumptions made in writing the 
truncated gap equation that are at fault. One obvious possible explanation is the absence 
of vertex corrections, although there are other missing skeleton diagrams which should also 
be investigated. 

The gap A is of course a gauge invariant observable. A complete calculation would yield 
a gauge invariant expression for the function /, which could be expanded as a power series 
in g. We learn three things from our (incomplete and gauge dependent) calculation. First, 
our results obtained in different gauges appear to converge at small g and support previous 
estimates of \img^Qf{g), namely the (7° term in the expansion of /. Second, because the 
results we obtain in different gauges only begin to converge for g < gc ^ 0.8, we learn that 
contributions to our gauge dependent function / which are of order g^ and higher must have 
gauge dependent parts which are numerically large at (7 ~ gc- Although we have simply 
evaluated f{g) and not expanded it in g, we learn that such an expansion is uncontrolled for 
g > gc- This suggests that if we knew the complete, gauge invariant function /, the g^ and 
higher terms in that expansion would also become uncontrolled for g > gc- It may be that the 
vertex corrections are the dominant contribution to the missing physics which is responsible 
for this breakdown: this hypothesis is supported by the arguments of Ref. that these 
effects contribute to / at order g^. Regardless of whether the vertex corrections turn out 
to be the most important effect left out of the truncated gap equation (|^, our calculation 
demonstrates that some contribution which is formally subleading is in fact large enough 
to render the calculation uncontrolled at (7 ~ gc- The third thing we learn is that although 
present calculations do yield reasonable estimates of lim^^o fig)^ if one is interested in using 
these calculations to estimate the value of A to within a factor of two, this can only be done 
for g <t:gc^ 0.8. 

In the CFL phase, all eight gluons get a mass. This means that in the CFL phase there 
are no gapless fermionic excitations, and no massless gluonic excitations, and therefore no 
non-Abelian physics in the infrared to obstruct weak-coupling calculations. The lesson we 
have learned is that even though everything is in principle under control, present weak- 
coupling calculations break down for g > gc ^ 0.8, corresponding to n < Hc with Hc ~ 10^ 
MeV (or higher [^|). This break down occurs even though A <^ ^ a.t g gc- It should be 
noted that what breaks down is the weak-coupling calculation of the magnitude of the gap 
A. Estimates based on models normalized to give reasonable zero density phenomenology 
can still be used as a guide, albeit a qualitative one. Furthermore, regardless of the fact 
that a controlled calculation of A has not yet been done at < 10^ MeV, it is possible 
to construct a controlled effective field theory which describes the infrared physics of the 
CFL phase on length scales long compared to 1/A, since in such an effective theory A is 
simply a parameter determined by physics outside the effective theory. This infrared physics 
is dominated by the massless Abelian gauge bosons ||T0|j39[| , the Nambu-Goldstone boson 
arising from spontaneously broken U{1)b 10, and the pseudo-Nambu-Goldstone bosons 
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arising from spontaneously broken chiral symmetry which have small masses due to the 
nonvanishing quark masses 1 10, 40 -47 . 
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APPENDIX A: THE MEISSNER EFFECT 

In this appendix, we set up the calculation of the Meissner effect. That is, we investigate 
the effect of the presence of a gap A on the functions F and G which describe the screening 
of the gluon propagator. 

In order to establish some necessary notation, we must begin by filling in some details 
in the derivation of Eq. (P) from Eq. (|^. We work in a color- flavor basis ({i, a}, {j, h}). In 
this basis, we define the following two 9x9 matrices: 



/o 
1 
1 



-1 



-1 
-1 





-1 



-1 





(Al) 



ab 



/2 1 



V 



1 

1 



1 

0/ 



(A2) 



which represent the antisymmetric color and flavor 3a and the symmetric color and flavor 
65 channels respectively in this basis. 
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In the derivation of the gap equation, we were only interested in the off-diagonal lower 
left component of the Nambu-Gorkov fermion propagator S. However, the calculation of the 
Meissner effect involves all components of the fermion propagator. Obtaining the fermion 
propagator by inverting the inverse propagator (Q) is straightforward but tedious. After a 
lot of algebra and using the ansatz (H) for the gap matrix, we find: 



S{q) 



Su{q) Si2{q) 
S2i{q) S22{q) 



(A3) 



where 



/ Aiq) B{q) B{q) 
B{q) A{q) B{q) 
B{q) B{q) A{q) 



Sii{q) 



C{q) 



C{q) 



C{q) 



C{q) 



C{q) 



\ 

( E{q) H{q) H{q) 

H{q) E{q) H{q) 

H{q) H{q) E{q) 



(A4) 



C{q)) 



S22{q) 



D{q) 



D{q) 



D{q) 



D{q) 



D{q) 



V 

(K{q) L{q) L{q) 

L{q) K{q) L{q) 

L{q) L{q) K{q) 



(A5) 



D{q)J 



S2iiq) = Si2iq) 



M{q) 
M{q) 



M(g) 
M(q) 



M{q) 
M{q) 



(A6) 



and where the above functions are defined as follows: 



16 



A{q) =7° 



p / X go-M-|g1 gg-(|q1+/.)^-3(A^(9„))^-ll(Af(9„)f-10A^(go)Af(go) 

+ 5g-(l9l+M)2-4(A^(5o)+2A^'(go))2 ^^T^l+^^F^(AjtaO^AfOT 



9o'-(l9l-M)'-4(Af(go)+2Af(go))2 



9g-(|g1-M)'-3(Af(go))'-ll(Af(go))'-lOA^^(qo)Af(qo) 



'zg-(k1-M)2-(Af(go)-Af(Qo))2 



5(g) 



go-M-|g1 



(A^(9o)+5Af(9o))(A^(go)+Af(go)) 



<?o'-(l9l+M)2-4{A^(go)+2A2^'(</o))2 </o^-(|gl+M)2-(A^(9o)-A|(go))= 



, p / N_ 90-M+I9I (Af(qo)+5Af(qo))(Af(qo)+Af(go)) 

^ -^''^^-(k1-M)^-4(Af(qo)+2A«(qo))2 gg-(|9l-M)2-(Af (qo)-Af(go))2 



(A7) 



go-A^-|g1 



go-M+l'?l 

g^-(|gl-M)2-(Af(go)-Af(go))2 



^^^^o'-(l9l+/^)'-(A^(go)-A2^'(go))2 ^ ^+^'^^o'-(l9l-/')'-(Af(5o)-Af(go))2 



go+M-|g1 



C 



gu+;.'T|t7l 



+ ^+^^)g^-(|gl-M)2-4(Af(go)+2Af(go))2 



g(i-(l^-/0'--HA2^''/"))"-iH^2 (g<)))"-i"A2Hg<))A2''wo) 



^- 9g-(|gl+M)2-4(A^(qo)+2At(«o))2 gg-(lgl+/^)^-(A^(go)-Af(5o))2 

go+/.-|g1 gg-(|g1-M)"-3(A^^(go))^-ll(Af(g„))^-lOAf(go)Af(go) " 



gg-(|g1-/^)'-(Af(go)-Af(go))2 



g5-(lg1+M)^-4(A^(qo)+2A|(qo))2 gg-(k1+Ai)^-(A^Hg(.)-Af(qo))2 

qo+/.-|q-1 (Af(go)+5Af(go))(A^^(go)+Af(go)) - 

gg-(k1-/^)2-4(Af(go)+2Af(qo))2 gg-(|ql-M)2-(Af (go)-Af(go))2 



c 



p I Af(go) 

[ gg-(|gl+M)2-4(A^(go)+2A«(go))2 



A^(go)-Af(go) (A^(go)+5Af (gp)) (A,^(go)+Af (gp)) 

gg-(lg1+A')^-4(A^(go)+2A«(go))2 gg-(|g-l+,.)2-(A^(go)-A|(go))2 



> p (n-] f Af(gp) 

^ ^^^^ Vg^(l9l-M)2-4(Af(go)-2Af (go))2 



A^(gp)-Af(gp) (Af (gp)+5Af (gp)) (A^^(go)+Af (gp)) 

^gg-(lg1-/^)2-4(Af(gp)+2Af(go))2 g§-(|gl-/.)2-(Af (gp)-Af (gp))2 



p f Af(gp)+A^(gp) 

+ ^^^ l^gg-(|g1+/^)2-4(A|(go)+2A|(gp))2 

-A^(qo)+Af(go) (A^(go)+5Af(go))(A^(go)+Af(go)) \ 
^gg-(|gl+/^)2-4(A^(go)+2A«(9o))2 gg-(|gl+M)2-(A^(go)-Af (go))^ J 
(n\( Af(go)+A^^(go) 
^ ^^''^ Vgg-(l9l-M)2-4(Af{go)+2Af'(go))2 



-A^(gp)+Af(go) (Af (gp)+5Af (gp)) (A^ (gp)+Af (gp)) \ 

go'-(|gl-M)'-4(Af(go)+2Aj(go))2 gg^_(|gl_^)2_(Af (gp)-Af (gp))2 j 



M(g) =C7^ 



P I'n'l -A^(gp)+Af(gp) p / x -A^^(gp)+Af (gp) 

+ ^^^gg-(|g1+/^)'-(A^(9o)-A2*(go))2 ^ ^-^y^gg-(|gl-/.)2-(Af(go)-Af(go))2 



c 



(A7) 
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k 

FIG. 3. One-loop contribution to the Meissner effect. 



Note that 5'2i(g) = Suiq) is a general property of the Fermion propagator S and can be 
proved for an arbitrary number of colors and flavors using only the definition of the inverse 
Fermion propagator, Eq. (p, and properties of the Dirac gamma matrices. Whereas only 
K, L and M were used in the derivation of the gap equation, all these functions are required 
in evaluating the Meissner effect. 

The Meissner effect is the change in the screening of the gluon propagator induced by 
the presence of a gap. To one loop order, we need to evaluate the gluon propagator of Fig. 
^ using the full fermion propagator including the gap. The result can still be written in the 
form ([TI|) but now 



F(g) = Fo(g) + 5F{q) and = Go(g) + 5Giq) 



(A8) 



where Fq and Gq are the A = functions written as F and G in (|T^). Recall that Go, which 
describes Landau damping, vanishes for go ^ 0. Because 6G is nonzero in the go — > hmit, 
the Meissner effect can be described as giving a mass to the gluons. Previous analyses of 
the Meissner effect have either been done for two-flavor QCD or have used simplified 

Our goal is to formulate the correct calculation of 6F{q) and SG{q) in 



estimates 
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the CFL phase. Recent work along the same lines can be found in Ref. 
From the diagram of Fig. R, we obtain the gluon polarization 



V / #Tr 



d*k 



Tr 



r^Suik + q)Y^Suik) + [r^f S22{k + q) [^^f S22{k) 
-rkS,2{k + q) [Y^f S2i{k) - {r^fs2iik + q)Y^S,2{k) 



(A9) 



where the trace is taken over color, flavor, and Dirac indices and all four elements of the 
fermion propagator, S{q), have been defined previously in Eqs. ( [A4D - ( |A7| ). This polariza- 
tion amplitude contains all the one loop contributions to the gluon propagator including the 
gap independent contributions, Fo{q) and Go{q). Ilab can be written in terms of F and G 
in a simple fashion: 



Jab 



(Goiq) + 6Giq)) P^-"^ + (Fo(g) + 6F{q)) P'^^^ 



(AlO) 



Hence, we only need to compute two components of U'^^ in order to obtain the functions 
SF{q) and SG{q), for example, Il^g and UH- Because we already know Fo(q') and Go{q), 
our goal is to extract SF{q) and 6G{q). We are therefore only interested in the difference 
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n[^^(A 7^ 0) — n(^b (A = 0). Finally, because SF{q) and SG{q) depend only on go and we 
can choose q to lie along the z-axis for simplicity. Keeping all this in mind, we find that (in 
Euclidean space) 

5Fiq) = ^i±M! (noo(A ^ 0) - n°°(A = 0)) 

SGiq) = UlliAy^O)-Ull{A = 0). (All) 

Note that (unlike the integrals which arise on the right hand side of the gap equation) the 
integrals which must be done in evaluating n(g) are ultraviolet divergent, and therefore 
sensitive to how they are cutoff at large ko and k. This ultraviolet divergence has nothing to 
do with A, and is canceled in our calculation of 6F and 6G by subtracting the A = result 
for n(g). We have checked that our results for 6F and SO are insensitive to the ultraviolet 
cutoffs in the integrals. 

Looking back at the definition of n(^^, we can see that it depends on Af'^{ko) and 
A^''^(fco). We make the same assumptions here as in our solution of the gap equation, 
namely that the antiparticle and sextet contributions can be neglected if A ^ /i and if one 
is interested in physics dominated by particles and holes near the Fermi surface. Before we 
proceed, let us define the following notation for the functions A{q) through M{q) defined in 
Eq. (|A7|) : identify the scalar functions multiplying the P± projectors with the appropriate 
± signs, e.g. A+{q). With this notation, the dominant contributions to the two polarization 
amplitudes we are interested in are: 

= -i^V(0r (l + {kTq) ■ k) [A.{k + q)A^{k) - B^k + q)B.{k) 

+2C_{k + q)C.{k) + E+{k + q)E+{k) -H+{k + q)H+{k) 
+2D+{k + q)D+{k) - 2K_{k + q)K.{k) 
+2L^{k + q)L_{k) + 2M_{k + g)M_(A;)] 

ni = -|^?'/(gr (l + 2{k^qik^ - {k^q) ■ k) [A_{k + q)A_{k) ^^^'^^ 

-B_{k + q)B_{k) + 2C_{k + q)C_{k) + E+{k + q)E+{k) 
-H+{k + q)H+{k) + 2D+{k + q)D+{k) + 2K_{k + q)K_{k) 
-2L_{k + q)L^{k) - 2M_{k + q)M_{k)] . 

In any one gauge, i.e. for a particular choice of ^, our task is now clear. We first calculate 
A(fco) with SF{q) = SG{q) = 0, as described in the body of the paper. We must then use 
( |A12| ) to evaluate SF{q) and SG{q) given by Eq. ( |A11D . As in the calculation of A, we can do 
all angular integrals analytically and evaluate the double integral over feg and | k \ numerically. 
We must then re-evaluate A(A;o) with the new gluon propagator, modified by the addition 
of SF{q) and SG{q). We must then iterate this procedure, calculating SF{q) and SG{q) and 
then recalculating A(/co) repeatedly, until all results have converged. We have not carried 
this program to completion. However, preliminary numerical investigation suggests that, in 
agreement with arguments and estimates made by others |2^J2^J2^J2^J5^J5^ , the change in 
A arising from the inclusion of SF and 6G is small. In particular, it appears to be much 
smaller than the change in A which arises if one changes gauge from ^ = — lto^ = Oto 
^ = 1. Perhaps at some extremely small g, the infiuence of the Meissner effect on the gap 
could be larger than the influence of the neglected physics whose absence we diagnose via the 
gauge dependence of our results. At any g at which we have been able to obtain numerical 
results, however, the Meissner effect is insigniflcant relative to that which is missing. 
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